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PREFACE 


The  author  was  very  much  attracted  by  the  very  fine  and  intense  beam  of  the  laser  when  he 
was  working  on  a  laser  interferometer  in  collaboration  with  Prof  Otto  R.  Frisch  at  the  Cavendish 
Laboratory,  University  of  Cambridge,  for  a  bubble  chamber  particle  track  measuring  device,  the 
SWEEPNIK.  The  very  fine  beam  enabled  the  track  measurement  to  be  accurate  to  a  few  microns. 
Later,  high  power  lasers  were  used  for  welding  and  cutting,  and  in  the  study  of  thermal  nuclear 
fusion.  The  diffusion  of  beam  energy  plays  the  central  role  in  all  these  applications.  The  diffiision 
of  beam  energy  comes  into  the  picture  again  in  the  study  of  skin  bum  and  eye  damage  by  lasers. 
This  report  concerns  the  basic  questions  regarding  the  boundary  condition  of  energy  diflEusion 
which  seems  to  be  the  key  to  understand  the  great  variety  of  diffiision  phenomena.  A  great  deal 
of  time  has  been  spent  to  program  and  properly  mn  the  newly  installed  mainframe  computer  at 
Natick  for  scientific  computing.  The  work  was  funded  as  AH98CA0B00,  Thermal  Protection , 
covering  the  period  from  November  1994  to  October  1995.  Programs  written  previously  were 
most  helpful. 
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THE  MATHEMATICAL  FOUNDATION  AND  THE  BOUNDARY 
CONDITIONS  ENCOUNTERED  IN  LASER  ENERGY  DIFFUSION  IN 
A  SOLID  AND  SKIN  SIMULANT  TEMPERATURE  RISE 


1.  Introduction 

The  diffusion  of  energy  from  a  laser  beam  is  a  very  important  process  in  laser  physics 
and  relates  to  a  great  many  areas  of  application  of  a  laser,  from  the  protection  of  skin 
and  eye  of  personnel  against  such  radiation  and  the  many  applications  in  medical 
science  of  a  low  energy  laser,  to  laser  welding,  laser  cutting  and  the  like  by  a  high 
power  laser,  not  to  mention  the  role  the  energy  diffusion  plays  in  the  very  high  power 
laser-induced  thermal  nuclear  fusion. 

As  the  diffusion  of  beam  energy  relates  to  so  many  areas  of  laser  applications,  the 
problem  was  critically  analyzed,  and  the  solutions  were  calculated  for  a  skin  simulant 
exposed  to  laser  radiation  with  exposures  ranging  from  10  seconds  to  20  microseconds. 

As  will  be  seen  in  the  following  sections,  the  solutions  involve  two  sets  of  boundary 
conditions  which  are  contradicting  to  each  other.  This  then  leads  to  the  fundamental 
questions  regarding  the  boundary  conditions  of  laser  energy  diffusion,  the  thermody¬ 
namics  and  the  diffusion  equation  itself. 

2.  Thermal  Diffusion  and  Boundary  Conditions 

The  diffusion  of  heat  in  a  solid  is  given  by  the  thermal  diffusion  equation 

cp^-KV^  =  0  (1) 

where:  T  =  temperature  at  (x,y,z,t) 

t  =  time 

p  =  density  of  the  solid 

c  =  specific  heat  capacity  of  the  solid 

K  =  thermal  conductivity  of  the  soHd 
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V^=  Laplacian  operator 


In  a  number  of  cases,  one  can  simplify  the  diffusion  Eq(l)  to  diffusion  in  one 
dimension. 

Thus, 


dT  ^ 


(2) 


For  a  one-dimensional  irradiation  say,  with  a  laser  beam,  one  can  further  approximate 
the  situation  as  an  energy  stream  flowing  into  a  semi-infinite  sohd. 

Two  boundary  conditions  are  to  be  satisfied  in  such  a  case.  One  is  that  at  the  far 
end  of  the  solid,  there  will  be  no  temperature  rise,  i.e.. 


r  =  o 


when  X 


oo 


(3) 


Here  and  in  Eq  (2),  T  will  from  now  on  denote  the  temperature  rise,  i.e.,  the  initial 
temperature  of  the  solid  is  taken  to  be  constant  with  respect  to  time  t  and  position 
X.  The  other  boundary  condition  is  that  at  the  surface  of  irradiation.  One  way  of 
tackling  the  problem  is  simply  to  equate  the  energy  flux  density  of  the  irradiation  to 
the  expression  of  heat  flow  at  the  surface  inside  the  solid, 

^0  =  l„o  (4) 

at  a:  =  0. 

We  will  call  this  the  constant  surface  heat  flux  solution^"^. 

Theory  of  differential  equations  then  says  that  Eq  (2)  together  with  the  boundary 
conditions  (3)  and  (4)  will  determine  the  values  of  T  at  any  space-time  point  (x,t). 
The  physical  properties  of  the  material  of  the  solid  K,  c  and  p  enter  into  the  problem 
via  Eq  (2). 
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3.  Boundary  Conditions  for  an  Irradiating  Laser  Beam 

For  irradiation  with  a  laser  beam,  equating  the  energy  flux  of  the  beam  to  the  heat 
flow  is  not  very  convincing.  For  example,  if  the  solid  is  transparent  for  that  particular 
radiation,  there  is  then  no  energy  absorption  from  the  beam,  and  hence  there  will  be 
no  temperature  rise.  Thus,  absorption  has  to  be  taken  into  consideration. 

For  a  medium  with  an  absorption  coefficient  a,  the  beam  with  energy  flux  density 
$0  passing  through  it  will  be  attenuated  as: 


- =  — 

di 

(5) 

or 

$  = 

(6) 

for  the  path  length  1. 

The  energy  absorbed  will  eventually  turn  into  heat  in  the  medium 

diffusion  equation  ,  instead  of  Eq(l),  will  then  be: 

.  The  thermal 

cp^  =  KV^T+Q 

(7) 

where  Q  is  the  energy  absorbed  per  unit  volume  per  unit  time,  or  in 

the  one  dimen- 

sional  case, 

dT 

dx 

(8) 

or 

d^T  IdT  a^o 

K  dt  K  ^ 

(9) 

where  /c  =  ^  is  the  thermal  diffusivity  of  the  solid. 

The  boundary  condition  at  the  far  end  will  be  the  same  as  before, 

o 

II 

(10) 

when  X 


oo. 
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At  the  surface  of  irradiation,  there  is  no  heat  flow 


-K§:=o 

OX 


(11) 


at  a:  =  0. 

4.  The  Zero  Surface  Heat  Flux  Solution 

Eq  (11)  is  at  variance  with  the  boundary  condition  Eq  (4).  In  order  to  appreciate 
the  meaning  of  these  two  different  formalisms,  one  may  want  to  see  how  Eq  (8)  or 
(9)  is  arrived  at.  This  will  be  given  in  another  report  together  with  its  solution.  For 
now,  we  write  down  the  solution  of  Eq  (9)  with  boundary  conditions  (10)  and  (11) 
as  follows  ($0  is  the  energy  flux  density  of  the  laser  beam): 


^0 


1  $0 


a^Kt—cex 


2Ka' 


.erfc 


a 


— 


2y/^ 


where 


TOO 

ierfcy  =  /  erfctdt 

Jy 

and  erfct  is  the  complimentary  error  function 


(12) 


(13) 
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erfct  =  —j=  /  e““  du  =  1  -  trfi 

yJ'K  Jt 


(14) 


erft  =  [  e  du 

yJ'K  Jo 


(15) 


For  irradiation  of  duration  T,the  temperature  rise  for  0<  t  <  oo  will  be  given  by 


AT  =  U{t).T{x,t)-U{t-T).T{x,t-x) 


(16) 
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where  T(x,t)  is  the  function  T(x,t)  given  by  Eq.  (12),  and  U(t),U(t-T)  =Heaviside 
functions  at  t  and  t-r. 

In  the  original  solution^’®  (no  derivation  was  given)  of  Eq  (12)  a  factor  a.  in  the 
denominator  of  their  third  term  was  missing^.  We  verified  that  Eq  (12)  is  the  correct 
solution. 

5.  Computations  for  a  Skin  Simulant  -  Long  Duration  Pulse 

With  the  aid  of  a  computer,  the  solution  Eq  (12)  and  Eq.  (16),  can  be  computed  for 
different  cases.  For  skin  simulant,  the  results  obtained  may  be  characterized  by  two 
groups:  (1)  with  exposure  to  the  laser  from  a  few  seconds  to  say  10  milliseconds  and 
(2),  from  a  few  milliseconds  to  20  microseconds.  Figs.  1  and  2,  with  exposure  time  r  = 
10  seconds  and  0.5  second,  show  the  typical  heating  characteristics  of  the  exposures  to 
the  longer  duration  pulses.  The  curves  in  both  figures  are  the  temperature  rise  curves 
at  depth  x  =  0,40,80,120  and  160/um,  with  =  0.167  K  =  1.16x10  ^  cm^^c  ’ 
K  =  1.11x10“^^^,  and  a  =  840cm~^. 

One  notes  in  Fig.  2  the  peaks  occur  from  ;^  =  1.0  at  x  =  0,  to  ~  1.5  at  x  =  160/fm, 
and  the  pulse  shape  is  broadened  as  the  temperature  wave  penetrates  into  the  material 
(skin).  This  is  indeed  the  case  ,as  we  remember  that  the  temperature  wave  propagates 
with  a  velocity  depending  on  the  frequency  of  the  temperature  source  and  hence  one 
will  see  the  dispersion  or  the  broadening  of  a  driving  pulse. 

The  occurrence  of  the  temperature  peaks  at  various  depths  in  Fig.  2  is  shown  in 


Table  1. 


0.5  Sec.  Laser  Pulse  Exposure 
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Figure  2.  Temperature  rise  curves  of  skin  simulant  exposed  to  a  laser  pulse  of  x  =  0.5s,  <{>0=  0. 167  cal  cm 
at  depths  x  =  0,  40,  80,  120  and  160 /im 


Table  1 


TEMPERATURE  PEAKS  AT  VARIOUS  DEPTHS  OF  SKIN 
SIMULANT  EXPOSED  TO  LASER  PULSE  OF  r=  0.55  AND  $o  =  0.167 

CAL 

CM^S 


Depth  X  (microns) 

Time  of  Peak  - 

r 

0 

1.0 

40 

1.02 

80 

1.10 

120 

1.30 

160 

1.50 

6.  Millisecond  and  Microsecond  Pulses  -  Differences  in  the  Two 
Formulations  and  the  Depths  of  Penetration 

We  now  turn  to  irradiation  of  short  pulses. 

(a)  1  Millisecond  Pulse 

With  a  1  ms  pulse  and  energy  flux  density  =  40  cal  cm"2s~\  the  surface 
temperature  rise  at  the  end  of  the  pulse  is  18.8°C,  and  at  x  =  40^m,  the  temperature 
rise  is  l.T^C.  These  are  the  two  points  plotted  in  Fig.3  joined  by  the  heavy  solid  line 
(1).  With  $0  =  16.7  cal  cm~^s“^,  the  surface  temperature  rise  is  7.9°C,  and  at  depth 
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X  —  40/xm  ,  the  temperature  rise  is  0.7°C.  These  are  joined  by  the  heavy  dashed  line 
(2)  in  the  figure.  These  are  the  temperature  rise  values  calculated  from  Eq  (12),  i.e., 
from  the  heat  diffusion  equation  Eq.  (9)  with  zero  heat  flux  boundary  condition  Eq 
(11)  at  x=0.  The  corresponding  values  at  x=0  and  40/im  given  by  diffusion  equation 
Eq.  (2)  and  the  constant  heat  flux  boundary  condition  Eq.  (4)  at  x  =  0  are  41.0“C 
and  0.1°C  respectively,  for  $o  =  40  cal  cm“^s~^  (light  solid  lines  (1*)  in  Fig  3);  and 
17.1°C  and  0.05°C  respectively  for  $o  =  16.7  cal  cm“^5“^  (light  dashed  line  (2  )  in 
Fig.3). 

The  conventions  of  the  lines  plotted  here  and  in  the  following  figures  are:  Solutions 
with  zero  surface  heat  flux  boundary  condition  are  plotted  in  heavy  solid  lines  (labeled 
as  (1))  and  heavy  dashed  lines  (labeled  as  (2)),  solid  lines  for  40  flux  density, 
and  dashed  lines  for  16.7  Solutions  with  constant  surface  heat  flux  boundary 

conditions  are  plotted  in  light  solid  lines  (l’),  and  light  dashed  lines  (2*),  again,  sohd 
lines  for  40  flux  density,  and  dashed  lines  for  16.7 

It  is  to  be  noted  that  the  values  of  the  temperature  rise,  both  at  the  surface  and 
underneath,  given  by  the  conventional  constant  surface  heat  flux  formalism  of  Eqs 
(2),  (3)  and  (4),  are  higher  than  that  by  the  more  realistic  and  logical  consideration 
of  Eqs.  (9),  (10)  and  (11),  i.e.,  the  zero  surface  heat  flux  solution  (12).  In  fact,  as 
the  laser  pulses  become  shorter  and  shorter,  the  difference  becomes  larger  and  larger, 
and,  as  will  be  seen  below,  with  a  20/is  laser  pulse,  the  constant  heat  flux  formahsm 
gives  values  ten  times  higher  than  that  given  by  zero  surface  heat  flux  solution  (12). 
In  the  meantime,  solution  (12)  shows  that  the  laser  beam  penetrates  much  deeper 
than  that  predicted  by  Eqs  (2),  (3)  and  (4).  For  exposures  of  seconds  duration,  the 
difference  given  by  the  two  formalisms  is  small  and  may  be  neglected  for  all  practical 
purposes. 
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Temperature  Rise  in 


Depth  in  aM 


Figure  3.  Temperature  rise  of  skin  simulant  at  depths  x  =  0,  40,  and  80  /xmfor  laser 
exposure  of  t  =  1ms,  <t)o=  40  and  16.7  cal  cm'^s'^  Heavy  lines  are  with  3T/8x|;t^=0, 
light  lines  with  -K3T/5x|,(,o'=  4>o  boundary  conditions.  The  figure  shows  the  much 
higher  surface  temperature  rise  of  a  solid  given  by  the  constant  surface  heat  flux  solution 
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Figure  4.  Temperature  rise  of  a  skin  simulant  at  depths  x  =  0,  40,  and  80  (itti  for 
laser  exposure  of  t  =  400  ps,  4)o=  40  and  16.7  cal  cm'^s'*.  The  figure  shows  the 
much  higher  surface  temperature  rise  of  a  solid  given  by  the  constant  surface  heat 
flux  solution. 
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Figure  5.  Temperature  rise  of  a  skin  simulant  at  depths  x  =  0.  40,  and  80  fim  for  laser 
exposure  of  t  =  100  ps,  4>o=  40  and  16.7  cal  cm’^s’'.  The  figure  shows  the  much  higher 
surface  temperature  rise  of  a  solid  given  by  the  constant  surface  heat  flux  solution. 
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Temperature  Rise  in 


50  jis  Laser  Pulse  Exposure 

00  =  40  cal  cm  *2  s  *' 

00  =  16.7  cal  cm  -^s  •’ 

Heavy  lines  (1)  &  (2):  Boundary  condition  O 

Light  lines  (1')  &  (2’):  Boundary  condition  |i|  =  0o/(-K) 


Depth  in  pM 

Figure  6.  Temperature  rise  of  a  skin  simulant  at  depths  x  =  0,  40,  and  80  /imfor 
laser  exposure  of  x  =  50  ps,  16.7  cal  cm  ^  s  .  The  figure  shows  the 

much  higher  surface  temperature  rise  of  a  solid  given  by  the  constant  surface  heat 
flux  solution. 


Temperature  Rise  in 


Figure  7.  Temperature  rise  of  a  skin  simulant  at  depths  x  =  0,  40,  and  SOfim 
for  laser  exposure  of  t  =  20  ps,  4^0=  40  and  16.7  cal  cm’^s'*.  The  figure  shows 
the  much  higher  surface  temperature  rise  of  a  solid  given  by  the  constant  surface 
heat  flux  solution. 
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(b)  400/^5  Pulse 

For  400  fis  pulse  of  the  same  energy  flux  density  #o  =  40  cal  ,  the  tem¬ 

perature  rise  calculated  from  Eq  (12)  at  the  surface  is  8.9"C;  at  x  =  AOfxm  ,  0.5°C; 
and  at  X  =  S0/J,m,0.O2°C  (joined  by  the  heavy  solid  lines  in  Fig  4).  For  $o  =  16.7cal 
cm“^s“^at  depths  x=0,  40/xm  and  80^m,  T  is  3.7°C,  0.2°C  and  0.009"C  respectively 
(heavy  dashed  lines  (2)  in  Fig  4). 

Calculated  from  diffusion  equation  (2)  and  boundary  condition  (4)  at  x=0,  the 
corresponding  values  for  $o  =  40  cal  cm“^s“^  are  25.4°C  at  x=0  and  1.5°xlO”‘*C  at 
X  =  40)um  (joined  by  light  solid  lines  (P  ),  Fig  4);  and  for  $o  =  16.7  cal  cm“^s“^, 
10.8‘’C  and  6.4°xlO~®C  at  x=0  and  40/im  respectively  (joined  by  light  dotted  line 
(2' )  in  Fig  4).  At  80  /im,  the  temperature  rises  are  zero  (~10“^’^C)  for  both  energy 
flux  densities  and  are  not  plotted  in  the  Figure. 

(c)  100/zs,  50fis  and  20 fis  Pulses 

The  results  given  by  the  two  formalisms  for  these  short  fis  pulses  and  the  same 
energy  flux  densities  =  40cal  cm~^s“^and  16.7  cal  cm“^s~^,  are  shown  in  Tables 
2  and  3  along  with  those  of  1  ms  and  400  fis  exposures.  They  are  similarly  plotted 
in  Figs  5,  6  and  7. 

7.  Discussion 

As  remarked  above,  the  constant  surface  heat  flux  solution  gives  much  too  high  a 
temperature  rise  at  and  near  the  surface:  For  a  ms  pulse,  it  is  higher  by  a  factor  of 
41.0/18.9  or  «2.2;  for  a  400/is  pulse,  by  a  factor  of  2.8;  for  a  100/is  pulse,  by  a  factor 
of  4.9;  for  a  50//s  pulse,  by  a  factor  of  6.6;  and  for  a  20fis  pulse,  by  a  factor  of  9.8, 
than  that  given  by  the  zero  surface  heat  flux  solution.  On  the  other  hand,  the  depth 
of  penetration  of  the  laser  pulse  by  the  constant  surface  heat  flux  solution  becomes 
smaller  and  smaller  as  the  pulse  width  decreases.  This  behavior  as  depicted  by  the 
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TABLE  2 

TEMPERATURE  RISE  BY  THE  TWO  DIFFERENT  BOUNDARY 

CONDITIONS*  $0  =40 


T 

X  in  fiM 

$0=40  cal  cm  ‘^s  ^ 

O 

II 

o 

(f  |o=  -^o/K  ) 

1  ms 

0 

18.9C 

41.0C 

40 

1.69 

0.12 

80 

0.06 

7.2x10-^ 

400/is 

0 

8.93 

25.4 

40 

0.53 

1.5x10-^ 

80 

0.02 

0 

100/is 

0 

2.64 

13.0 

40 

0.12 

0 

80 

0.004 

0 

50fi$ 

0 

1.39 

9.17 

40 

0.06 

0 

80 

0.002 

0 

20^3 

iniiiii^^ 

0.59 

5.80 

IIIIQIIIIII 

0.02 

0 

80 

8.0x10-^ 

0 

*  Column  3  is  plotted  as  (1)  in  Figs  3-7 

*  Column  4  is  plotted  as  (l')  in  Figs  3-7 
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TABLE  3 


TEMPERATURE  RISE  BY  THE  TWO  DIFFERENT  BOUNDARY 

CONDITIONS*  $0  =  16.7 


*  Column  3  is  plotted  as  (2)  in  Figs  3-7 

*  Column  4  is  plotted  as  (2*)  in  Figs  3-7 

constant  surface  heat  flux  solution  is  not  surprising.  In  fact,  the  constant  surface  heat 
flux  boundary  condition  holds  only  in  the  ideal  case  in  which  the  surface  of  the  solid 


is  in  contact  with  a  constant  temperature  heat  reservoir  whose  temperature  follows 
in  every  small  (infinitesimal)  step  the  surface  temperature  rise  of  the  solid  or,  in  the 
case  of  irradiation,  when  the  beam  energy  (laser)  is  totally  converted  to  heat  at  the 
surface®’^.  The  penetration  of  heat  in  the  latter  case  is  due  to  diffusion  given  by  Eq. 
2  and  controlled  by  the  diffusion  coefficient  /c  =  ^.  It  is  to  be  noted  that  although 
in  the  figures  (Figs  3-7),  the  temperature  rise  versus  depth  graphs  are  all  drawn  as 
straight  lines,  they  are  actually  concave  . 

8.  CONCLUSION 

The  investigation  of  energy  diffusion  from  a  laser  beam  in  a  solid  has  produced 
many  useful  results  related  to  personnel  protection  against  such  radiation  and  many 
other  applications  of  the  laser.  It  was  presented  also  with  the  fiirida.Tnpnt.a1  questions 
regarding  the  boundary  conditions  of  the  diffusion  equations,  the  zero  surface  heat 
flux  vs  the  constant  surface  heat  flux,  ||  |j:=o=  0  vs  -K|^|j,=o  =  ^o,  at  the  boundary. 
The  constant  surface  heat  flux  boundary  condition  in  the  conventional  solution  of  the 
heat  diffusion  equation  is  shown  to  be  unrealistic  for  a  laser  beam.  Justification  for 
zero  surface  heat  flux  boundary  condition  for  the  laser  irradiation  is  presented  along 
with  the  calculations  for  the  temperature  rise  of  a  skin  simulant  exposed  to  the  laser 
irradiation.  The  former,  (constant  surface  heat  flux)  solution  gives  the  temperature 
rise  higher  by  a  factor  of  2.2,  2.8,  4.9,  6.6  and  9.8,  for  laser  exposure  of  1000,  400, 
100,  50  and  20  pts,  respectively,  than  that  given  by  the  diffusion  equation  with  zero 
surface  heat  flux  boundary  condition.  The  analysis  holds  for  the  general  case  of  other 
laser  irradiated  solids. 
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